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Pion exchange in S-wave between hadrons that are themselves in a relative S-wave is shown to 
shift energies by hundreds of MeV, leading to deeply bound quasi-molecular states. In the case of 
charmed mesons D* , Di a spectroscopy arises consistent with enigmatic charmonium states observed 
above 4 GeV in e + e~ annihilation. A possible explanation of Y (4260) — > iprnv and Y (4360) — > vp'-Kir 
is found. We give results for all isospin and charge-conjugation combinations, and comment on 
flavor exotic doubly charmed states and bottomonium analogs. A search in DDStt is recommended 
to test this hypothesis. An exotic 1 is predicted to occur in the vicinity of the Y(4260). 

I. INTRODUCTION 

If two hadrons A, B are linked by A — > Btt, then necessarily hadronic pairs of AB or AB have the potential to feel 
a force from n exchange. This force will be attractive in at least some channels. Long ago [HQ the idea of tt exchange 
between flavored mesons, in particular charm, was suggested as a source of potential "deusons" [l| . Using the deuteron 
binding as normalization, the attractive force between the J p = CP charmed D and its J p = 1~ counterpart D* was 
■ calculated for the DD +c.c. S-wave combination with total J = 1 ++ , and the results compared with the enigmatic 
charmonium state X(3872)0,|H|. 

In all these examples, as in the more traditional case of the nucleon, where the NNir coupling is the source of 
an attractive force that helps to form the deuteron, the exchanged tt was emitted and absorbed in a relative P-wave 
with respect to the hadrons. In such cases, the binding energies that result are 0(1 — 10)MeV; this in particular has 
encouraged interest in the X(3872) which is within errors degenerate with the D°D*° threshold. It has recently been 
pointed outQ that the exchange of a n in S-wave, between pairs of hadrons that are themselves in relative S-wave, 
leads to deeply bound states between those hadrons. Instead of binding of a few MeV, as in the cases considered 
historically, there is now the potential for binding on the scale of O(100)MeV, leading to a rich spectroscopy of states 
that are far from the di-hadron channels that create them. We shall argue that examples of such a spectroscopy 
& , appear to be manifested among charmonium-like mesons. 

We organize this paper as follows. First we summarize the general arguments for expecting large binding energies 
due to S-wave 7r exchange. We shall consider a chiral Lagrangian model to illustrate and quantify the phenomenon 
of energy shifts of O(100)MeV, focusing specifically on the charmonium-like 1 isoscalar (7 = 0) channel. In 
Section IIIII we investigate the connection between the chiral potential and the decay width of the relevant charmed 
mesons, first in the heavy quark limit with point particles and subsequently in the non-heavy quark limit and with 
form-factors. Then, we solve the Schrodinger equation and discuss the uncertainties within the model. Detailed 
results for the charmonium-like = are presented in Section llVl along with results for other J PC , isospin, 

and flavor channels. We discuss the limitations of our approximation to the strong interaction in Section |Vl give 
phcnomenological implications and suggest experimental searches in Section IVT1 and finish with conclusions in Section 
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II. MOLECULES AND S-WAVE tt EXCHANGE 



Several groups have studied the following meson-pairs looking for bound states in total J PC channels due to pion 
exchange (from here on DD etc. will be taken to include the charge conjugate channel): 

D*(l~) -> D(Q-)ir leading to the deuson DD* A(3872) J PC = 1++ 
D*(l~) -> D*(1~)tt leading to the deusons D*D* J PC = 0++, l+-,2++ 

These combinations were discussed in In all of these examples parity conservation requires that the tt is emitted 
in a P-wave; the hadrons involved at the emission vertices have their constituents in a relative s-wave, (we use S,P to 
denote the angular momentum between hadrons, and s,p to denote internal angular momentum of the constituents 
within a hadron). In such cases, the tt emission being in P-wave causes a penalty for small momentum transfer, q, 
which is manifested by the interaction Q 



q r + M 



where ^ 2 = m 2 — (tob — m^) 2 , being the masses of the mesons in the process A — > i?7r. (For a discussion of 

this interaction, and its sign, see Eq. (20) in refQ). The resulting potential is oc q 2 for low momentum transfer and 
has been found to give bindings on the scale of a few MeV, which is in part a reflection of the P-wave penalty. 

There is no such penalty when tt emission is in S-wave, which is allowed when the hadrons A, B have opposite 
parities. Examples involving the lightest charmed mesons are Di{l + ) — > D*(1~)tt and Z?o(0 + ) — > D{Q~)tt. P-wave 
tt exchange carries a q penalty at each ttAB vertex. One might naively anticipate that the transition from a D or 
D* , with constituents in s-wave, to the D\ or Dq, with constituents in p-wave, would restore a q penalty, leading to 
small binding effects as in the cases previously considered. However, as we now argue, this need not be the case, and 
energy shifts of O(100)MeV can arise. 

There is a long history of data on tt transitions between hadrons of opposite parity which indicate that the S-wave 
coupling is significant when q — >• 0. In the charm sector of interest here, the large widths^ for r(L> ->■ Dtt) ~ 
260±50MeV and r(L>i(2430) -> D*tt) ~ 385±O(100)MeV suggest that, even after phase space is taken into account, 
there is a significant transition amplitude. This non-suppression was specifically commented upon in the classic quark 
model paper of ref. [§[ . It arises from a derivative operator acting on the internal hadron wave function, which enables 
an internal s to p transition to occur without suppression even when q vanishes. This can be seen when ipjsip is 
expanded to <r.(q — wp/m), where p the internal quark momentum^. Feynman, et aZ.Q argued for this form on 
general grounds of Galilean invariance. The presence of p gives the required derivative operator, and hence the 
unsuppressed p s transition. 

An unsuppressed transition, when applied to tt exchange in the D\D* system (e.g. [Ioj |) causes the chiral model 
analogue of Eq.([T]) in the / = 0, channel to become 



where hj (y/2f n ) is the D\D*tt coupling constant (up to a phase), = 132 MeV, q is the exchanged three-momentum, 
fi 2 = —(mjj 1 — mjj*) 2 + to 2 (/1 2 < for the D\D* system), and (tj • Tj) is the usual contraction of Pauli matrices 
resulting from the exchange of an isovector by two isospin-half particles. T is the model dependent form-factor which 
regulates the potential and would be unity in the chiral model. 

In the derivation of the potential a static approximation has been made to the pion propagator. The full propagator 
is q 2 — m 2 = (Ea — Eb) 2 q 2 rn\. Approximating Ea — tua and Eb — tub one recovers the form of the propagator 
presented in the potential, Eq. ([2]). 

The potential is similiar to one presented in Table 1 of Ref. 10], who were investigating a tt exchange model of the 
1~, / = 1 Z + (4430). They considered only the 7 = 1 channel and omitted the Tj • Tj factor (which is unity for I = 1). 
We have made this factor explicit as it will become crucial when we study the I — sector later. The absence of a q 2 
penalty factor is immediately apparent. The scale is now being set by the mass gap squared, which is equal to the 
timelike component of the momentum transfer vector squared, (?g — > (mjj 1 — mjj*) 2 as |q| — > 0. 
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A -» Bn m A /MeV m a /MeV F/MeV BF |q|/MeV 

Di(2430) -> £>*(2010)tt 2427 ± 40 2010.27 ± .17 384 +^° N/A 359 

L»S(2400)° ->• Dtt 2352 ± 50 1896.62 ± .20 261 ± 50 N/A 414 

#5(2400)* ->■ Dn 2403 ± 40 1864.84 ± .17 283 ± 40 N/A 461 

Si (5721) -> B*(5325)vr 5723.4 ± 2.0 5325.1 ± .5 N/A dominant 360 

#1(1400)* -»■ if*(892)7r 1403 ± 7 891.66 ± .26 174 ± 13 94 ± 6% 402 

Ko (lMQ)* -> Kir 1425 ± 50 493.677 ± .016 270 ± 80 93 ± 10% 619 

TABLE I: Data of low lying mesons of different flavor sectors with opposite parity which exhibt a large width. Values taken 
from the Particle Data Group[3J. No width data are available for the bottom sector and no branching fractions are given for 
the charmed sector. 



This potential and any bound states have immediate implications for a rather rich set of physics. The potential also 
applies for the DqD system, and the bottomonium and strangeonium analogs of D\D* and D$D, by exchanging the 
masses with their appropriate counterparts. Note that the potential in Eq. ([2]) has no spin dependence and therefore 
any results apply equally to the D\D* spins coupled to total spin 0, 1 or 2. For example, if an isoscalar 1" bound 
state is found, then we also expect degenerate and 2~~ states. 

Thus on rather general grounds we may anticipate significant energy shifts, ~ O(100MeV), due to it exchange at 
least in some channels between such hadrons in a relative S-wave. Signals may be anticipated below or near threshold 
in the following channels (in the charmonium analogues, involving charm and anti-charm mesons for either / = or 
1, or in exotic states with manifest charm involving two charm mesons): 

L>o(0 + ) -> L>((T)7r leading to the deusons DD Q J PC = _± 
L>i(l+) -> L>*(1")tt leading to the deusons D*~~D 1 J PC = (0, 1,2)-* 



We also find that it is possible that L > states could bind which would lead to more J p configurations. 

Pion exchange depends on the presence of u, d flavors, therefore there will be no such effects in the D S D S analogues. 
Further, the potential depends only on the quantum numbers of the light quarks. Therefore, there will be effects in 
the strange and bottom analogues, which can add to the test of such dynamics at different kinematics. 

The parameter h in Eq. $2$ is closely connected to the width of the D\ — > D*ir decay. Data exists on this decay 
which constrains the value of h and hence the spectrum of the model. We discuss the extraction of h from the decay 
width in the next section. 



III. THE COUPLING CONSTANT h 



Being simply related to the DiD*ir coupling constant, h also appears in the chiral formula for the D\ — > D* decay 
width. Eq. (137) of Ref. [H gives: 



T(D» -> D*+n-) = i- (^) 2 - ™ D *f. (3) 

which is valid in the heavy quark limit. 

In order to extract h using Eq. Q , we use the data from the PDG listed in Table |U In the absence of a branching 
fraction we assume that the total width is saturated by the D*w channels. We are using chiral formulae for the charged 
7T width which may be related to the total tt decay width by T(Df -> £>*+7T~) = |r(£»? -> D*tt) [llj. Therefore, we 
use = 140.MeV and U = 132MeV throughout. For the D x -> D*ir system we have h = 0.63l;i|. 

There are theoretical and empirical reasons to suspect that Eq. ^ may be a poor estimate for h given T. Firstly, 
in the heavy quark limit assumed for Eq. (JH]), muj = m£> , mw = mp, and thus T(Di D*ir) — T(Do — > Dw) as 
Eq. applies equally well to the Dq — > Dir decay. However, these relations do not hold experimentally. Finite mass 
effects (including mass differences) have been used to derive a correction to Eq. ©[Hj]: 

+ _ h 2 \q\m D * 2 2 2 1 / (m Dl + m D .) 



8tt/2 m 3 D u 3 V Am 2 D m 2 D , » 
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System HQ NHQ NHQFF 

Di(2430) -» D*(2010)tt 0.63t;^ 0.80l 2 ? 7 l-Otaa 

£>S(2400)° Dvr 0.41 ± 0.06 0.55 ± 0.08 0.79 ±0.11 

7^2(2400)* ->• D-k 0.36 ±0.04 0.50 ± 0.05 0.80 ± 0.08 

Ki ( 1400) ± ->• iT (892) 7T 0.30 ±0.02 0.50 ±0.04 0.70 ± 0.05 
7^(1430)* -> ^tt 0.15 ±0.04 0.47 ±0.13 1.2 ± 0.3 

TABLE II: Values of h extracted for various systems which may experience S-wave n exchange. The adaptation of equations 
for the charm-system to their appropraite form for flavor analog systems by making obvious mass substitutions is assumed. 



Using this expression we have h = 0.80* q® 7 . 

We have mentioned that our analysis of the D\D* system applies equally well to the DqD system. Indeed, Eq. ^ 
applies to both systems with a trivial substition of the appropriate mass. However, when finite mass effects are 
included, chiral model gives a different formula for the decay widths of the Di and D mesons. The analogous 
formula to Eq. Q is[l2| 

rvn , n+ — \ h2 \<l\ m D , 2 2 N2 / k s 

T(A) -> D+tt ) = (m Do - m D ) . (5) 

Due to the larger mass gap (and hence the larger |q|), Eqs. (@| and ([5]) imply that T(D — » £>7r) 1.5r(Z3! — > D*tt). 
Empirically Q, 

T(D -> L>tt) - 260 ± 50 MeV and r(£>i (2430) -> L>* tt) - 385 ± 0(100) MeV. 

Although the uncertainties are large, T(Dq — » Dir) has a smaller width than T(Di — >• Z?*7r) even though the phase 
space is larger, in contrast with the expectations of Eqs. (j4]) and ([5]). 

In general, processes such as T(D® — » D* + n~) involve form factors that summarize the penalty for selecting the 
exclusive process of single it emission, which is increasingly improbable at large |q| relative to multi-pion, inclusive, 
channels. Thus, the assumption that the DiD*tt (DqDt:) coupling is constant in the chiral model does not take 
account of the full dynamics at the vertex. The data suggest that we must include the effects of exclusive form 
factors. 

The effects of form factors may be modelled by making the replacement h —> hJ-(\q\) everywhere. .F(|q|) is a 
smooth, decreasing function such that .F(|q| = 0) = 1. The exact form of T is model dependent; however, the 
introduction of a form factor will in general lead to an increased estimate of h and so, naively, to an increased binding 
energy. 

As an explicit example, consider the form factor resulting from a dynamical model of it emission [13j: 

F{x) = (l - e- 2 / 12 (6) 

with x = |q|//3 and /3 - 0.4GeV[13j. For D\ -> D*+ir- one has x = 0.89 while for D% -> D+tt" x = 1.18. This plays 
a significant role in the relative widths as 

^ = 0.89) l2 = L6 (?) 



LT{x = 1.18). 
which drives the ratio of widths in favour of the D\ . 

In turn the form-factor also shows that h, extracted earlier from the chiral model, is an underestimation. In such 
a model the more general Eq. (|4|) modified the heavy quark value of h = 0.63l; 13 to h — 0.80t; 2 ° and the effect of 
form- factors further increases h to h = l.Olg '\. Therefore, the inclusion of finite mass corrections and the effects 
of exclusive form factors has a significant impact on the value of h extracted from experimental decay widths. We 
emphasise that although the form factor itself is model dependent, the suppression for larger |q| is expected in general. 

In summary, from these different determinations we find values ranging from h ~ 0.5 to 1.3: the value of h is highly 
model and data dependent. We collect these results and present other results for analogous systems in Table HU The 
HQ column presents the values of h extracted in the heavy quark limit using Eq. ©. The NHQ column is similiar 
but extracts the values in the non- heavy quark limit using Eqs. (|4]) and The NHQFF column presents extracted 
h values which would be required to overcome the form- factor suppresion, Eq. ^ , and to reproduce the correct width 
in the non- heavy quark limit. In the following section we will present results for a range of h and show that the 
spectrum is highly sensitive to the value of h. 
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IV. MOLECULAR SPECTROSCOPY 



Previously @ we performed a variational calculation with the potential in Eq. (J2J and h = 0.8 ±0.1 using trial wave 
functions. With this technique we agreed with Ref. [l(| that there was no reason to expect an isovector \~~ ,D\D 
bound state. Additionally, we found deep binding in the isoscalar D\D system. The presence of deep binding in 
the 1 channel motivated the present study where we solve the Schrodinger equation and analyze the spectroscopy 
emerging from S-wave 7r exchange binding of the D\D and analogous systems. 

We solve the Schrodinger equation and quantify the bound states from S-wave ir exchange using a range of h to set 
the scale. The resulting spectrum contains several potential bound states. The 1 , / = channel includes IS and 2S 
states which are consistent with the V(4260) and V(4360) structures claimed in e + e~ annihilation. Results for the 
charmonium-like exotic 1 and isovector 1 channels are also presented. We find the binding energies are highly 
sensitive both to the value used for h in all channels, and attempts to model finite size effects in some channels. We 
first consider the potential from the chiral model involving a pointlike interaction, and then discuss modification of 
the potential due to finite size effects. 



A. Point-like pion exchange 



The Fourier transform of Eq.@ gives the D\D 1 potential with S-wave 7r exchange in coordinate space. When 



-jl 2 < the real part of the potential is: 



V s (r) 



h 2 {mrt 1 — mo*) 2 cos(pr) 



(8) 



in agreement with Ref. [T(J. We numerically solve the Schrodinger equation using this position space potential as 
described in Ref. Q. 

Clearly the results for larger values of h, which yield significant binding, will have important finite size corrections. 
Therefore the point particle results can only be considered to give a cursory quantitative examination of the impli- 
cations of our general argument. However, given the unusual nature of the oscillatory potential, it is beneficial to 
study the unregulated potential in order to contextualize the effects of a form-factor which are explored in the next 
subsection. 

In Table Hill we show the binding energies of some of the low-lying isoscalar D\D* states in relative S-wave with 
C = — (the parity obviously depends on the relative orbital angular momentum of the system; since the potential is 
independent of spin, the binding energies are degenerate across all possible total J combinations). Binding energies 
are given for a few values of h G [0.8, 1.3]. The binding energies are seen to be very sensitive to the value of h. 



Binding Energy / MeV 



State 


h = 0.8 


h = 0.9 


h = 1.0 


h = 1.1 


h = 1.2 


h = 1.3 


IS (0,1,2) — 


230 


415 


680 


1000 


1500 


2100 


2S 


12 


20 


29 


39 


76 


210 


3S 


1.5 


3.6 


6.7 


11 


51 


65 



TABLE III: Binding energies for various isoscalar D\D states in L 
values of h in the range identified above. 



with C — — ; the binding energies are given for a few 



If for example h = 1.0, we find that two or even three S-states may arise, with binding energies 680 MeV (IS), 
29 MeV (2S) and 7 MeV (3S). The rms radii, r rms for these states are then approximately 0.2 fm, 3 fm and 7 fm 
respectively. This shows that the ground state is typically hadronic, the 2S consistent with a canonical molecule and 
the 3S dubious. 

Using Fig. [TJ we can interpret the results for the r rms values obtained for the S-wave states with the point particle 
potential. The IS state had an r rms w 0.2fm clearly indicating that the state is bound in the first attractive well. In 
contrast, the 2S state had an r rms « 3fm indicating that the particles are bound in the second attractive well of the 
potential. The 3S state has an r rms ss 7fm suggesting that it is bound by the third attractive well. 

If we take the potential, Eq. and applied it to L > systems unchanged apart from the centrifugal potential, 
we would find multiple bound states in the P- and D- waves including some with binding energies of O(50)MeV. Firm 
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conclusions regarding the possiblity of such states would require a more extensive analysis of the origin of Eq. ([2]) than 
presented here. 

The potential energy scales ~ h 2 but the binding energies are much more sensitive to h (ground state binding 
energies scale like ~ h 6 ). This sensitivity to h may not be unexpected, as the oscillatory nature of the potential in 
position space makes the potential turn over to repulsive when r > 0.7fm, and gives considerable sensitivity to these 
oscillations even for the short range IS level, and critically so for the 2S. In a Coulomb potential the binding energies 
would scale as ~ h 4 ; this further explains the sensitivity noted above in the numerical calculation. 



B. Form Factor 



As noted previously, the form factor has a significant impact on the calculation of h from the decay width. In the 
previous section we used this "form-factor-renormalised" value of h, but otherwise continued to treat the potential as 
if the hadrons were pointlike. Therefore, it is prudent to investigate what effect form factors may have on the analysis 
of the molecular spectroscopy. To examine this question we attach dipole form factors to the potential, Eq. ([2]). Such 
ideas have been discussed in refs[l| ,[lH and (l6j . Following those ideas, we specifically choose to parametrise the form 
factor as 



T = 



A 2 - 



A 2 



A 2 



A 2 +H 2 



(9) 



and the potential is multiplied by T 2 - we use the latter expression for T . We have made the same static approximation 
as in Eq. ([2]). In nuclear physics dipole form factors have only A 2 + q 2 in the denominator due to the 7r exchange 
between the nearly degenerate nucleons. 

In position space, the form factor changes the potential from that in Eq. © to: 



V s {r) 



h 2 (mp 1 - mp*Y 



cos(/2r) 



,-JCr 



(A 2 



2X 



m l) e -Xr 



(10) 



with X 2 = A 2 + /i 2 — 771^2 = A 2 — (niDj — TO£>.) 2 and where /i 2 = — jl 2 < 0. We plot this potential for the isoscalar 
1 — - channel for various A in Fig. [1] and for a fixed A and the various l - channels in Fig. [SJ 

Here A is a purely phenomenological constant. Although its value should be related to the convolution of the 
spatial wave functions of the hadrons, its value is fairly arbitrary in practice. In the data-rich nucleon-nucleon 
sector, dipole form factors have been used in the Bonn nucleon-nucleon potentials. In CD-Bonn one finds values of 
A = 1.3 — 1.7GeV[l4|. However, there is no reason to believe that the value used in nuclear forces should be related 
to the value most appropriate for use in it exchange between charmed mesons. In the literature, other practitioners 
using dipole form factors in meson exchange molecular models employ values of: «1.2 GeV[l], ~ 1.2-2.3GeV 15], and 
« O.4-lOGeV0. 

The qualitative effects of this form-factor are made apparent in Fig. [1] Regulating the potential introduces a soft- 
repulsive core instead of a singular attraction at the origin. As A decreases, the first attractive well in the potential 
is entirely overwritten as a repulsive core. 

We present the results for the binding energy as a function of h and A for the IS and 2S isoscalar 1" D\D* states 
in Figs. [2] and [3] The horizontal axes are 1/A so that the origin corresponds to the point-like case. As one can see, 
the ground state binding energy falls off rapidly with decreasing A. Eventually the ground state binding energy finds 
a stable point and remains at approximately that energy for the rest of the considered values of A. This behavior is 
sharply contrasted by the binding energy of the 2S state. The 2S binding energy is initially insensitive to a decrease 
in A before falling steeply and finally becoming insensitive again. 

This behavior can be understood from the behavior of the potential in Fig.[TJ As A is decreased from oo, the potential 
is increasingly regulated. This manifests as overwriting the initial attraction from the potential and eventually 
replacing it by an entirely repulsive core interaction for A < 800MeV. Thus we would expect a steep fall off in ground 
state binding energy as A is decreased. This expectation is borne out in Fig. [2] In contrast, the 2S state is bound 
primarily by the second attractive well, which is unaffected by decreasing A as long as A > 800MeV. Thus, we would 
expect the 2S binding energy to be relatively stable against decreasing A as Fig. [3] confirms. At some point, which 
is h dependent, there will no longer be enough attraction in the first attractive well to bind the system, and so the 
ground state will begin to require presence in the second attractive well in order to bind, displacing the 2S state and 
decreasing its binding energy. When the first attractive well is completely overwritten, both the IS and 2S states 
should have a relatively stable binding energy as the attractive wells (second and third) which bind them are stable 



200 




r/fm 



FIG. 1: The potential, Eq. (|10[) . plotted against r in the isoscalar 1~~ channel for h — 0.8 and a variety of As. The solid line is 
the point particle result - in effect A — > oo; the dashed line is the result for A = 1.5GeV; the dash-dot line is for A = 1.0 GeV; 
the dash-dot-dot line is for A = .75GeV; and the dotted line is for A = .5GeV. 



against decreased A. Indeed the binding energies of the IS state decrease slightly as A — > 500MeV corresponding with 
the alteration of the second attractive well in Fig. [TJ 

This analysis shows that the molecular spectroscopy is very sensitive to the parameters. While a simple abstraction 
of parameters from existing data support the idea that a spectroscopy of molecules could arise, one cannot with 
certainty predict this. However, the result of strong binding appears relatively robust. Indeed our results show that 
the existence of robust bound states (assuming h is sufficently large) does not depend on deep attraction at the origin, 
and that, even in the presence of a strong repulsive core interaction, a bound state should exist with a binding energy 
largely determined by long-range 2fm) virtual pion effects. 

If the y(4260) and V(4360) are confirmed as genuine signals, then within this simple modelling, their energies are 
qualitatively consistent with those expected for S-wave binding. Indeed, the differing sensitivities of the IS and 2S 
states to A would allow one to tune the model to reproduce the binding energy of the V(4260), 174 ± 9 MeV and the 
y(4360), 76 ± 13 MeV assuming the mass of the Di was exactly 2427 MeV. Since the mass of the D\ affects both 
binding energies in a systematic way, we cannot simply add its error in quadrature for both to obtain our binding 
energies with their error. Instead, we study the system for m,D 1 = 2427 — 40MeV; 2427MeV; 2427 + 40 MeV requiring 
binding energies of: 134 ± 9MeV, 36 ± 13MeV; 174 ± 9MeV, 76 ± 13MeV; 214 ± 9MeV, 116 ± 13MeV. We present the 
"tune-ability" of the model in Fig. HJ 

The mass of the D\ effects the potential in two straightforward ways. First, it factors into the calculation of h 
from the decay widths. Secondly, it helps determine the mass gap which, along with h controls the strength of the 
potential. Although the value of h and the mass gap depend on the value of the D\ mass, 

hvs = h—— (11 

2427 — m o* 

is unchanged as varies over its error. This allows us to plot the different mass cases on a single axis. (The mass 
of the D* has an insignificant error.) 

Fig. [4] was produced by parameterizing the binding energies. We assumed that the 2S binding energy was approx- 
imately independent of A and so could be used to determine h. This assumption has been explicitly verified for the 
values of h, A considered and is found to be a good approximation. Then the IS binding energies were parameterized 
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FIG. 2: Plot of the IS 1 isoscalar binding energy for multiple values of h as the form factor parameter A is varied. 
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FIG. 3: Plot of the 2S 1 isoscalar binding energy for multiple values of h as the form factor parameter A is varied. 

as a quadratic function of 1/A whose coefficients were quadratic functions of h. This parameterization reproduced 
the IS binding energies over the relevant range of these parameters. The quadratic formula could then be used to 
extract the range of A from the F(4260) binding energy at each h. 

The region of compatability extends to just below the error bounds for /i e g to slightly above it. A values are 
undetermined by experiment, however the compatible values lie around 1 GeV which is near values used by other 
practictioners. Thus, a consistent, physically reasonable parameterization of h and A is possible which permits the 
identification of the Y (4260) and y(4360) as the IS and 2S bound states of the D\D system respectively. 

In general within the chiral model, ft, is a function of the experimental T(Di — > D*tt). If experiment were to show 
that the width were different than the current values that we have used, the consequent alteration of the molecular 
binding energies could be considerable. It is here that some major limitations in the robustness of the model lie. 
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FIG. 4: Countour plot of the values of h^s and A. The interior of the boxes corresponds to values which reproduce the binding 
energies of the y(4260) and F(4360) to within errors. The center box is for the experimental Di mass while the box on the 
left is for the Di mass minus its error and the box on the right is the D\ mass plus the error. The dotted line corresponds to 
the value for /ieh = 1.0 + 0.3 from the Di experimental width. 



C. Other Di D* bound states and flavour exotics 



The potentials in the C = ± and isovector/isoscalar channels are related by a simple constant. The potentials of 
isovector and isoscalar channels are related by a Tj • Tj factor while the potential in channels with opposite charge 
conjugation are related by a relative phase. Therefore, we can use Fig. [1] to intcrpet the binding in all these channels 
against finite size effects and a possible repulsive core. The robustness of our results in the isoscalar 1 channel were 
discussed previously. 

In Table ITVl we show the binding energies of some of the low- lying isoscalar and isovector D\D states in relative 
S-wave with C = ±. Binding energies are given for a few values of h € [0.8, 1.3] and A = 1 GeV. Interestingly we find 
potentially robust binding in all isospin and charge-conjugation states. 

We note that the pattern of binding described here is valid for A = 1 GeV and the pattern will be altered as A 
changes. In particular, the pattern will change if the finite size effects wipe away less of the deep attractive core which 
binds the isoscalar 1 and isovector 1 h channels. In general a higher value of A will lead to (significantly) more 
deeply bound isoscalar 1 and isovector 1 h bound states, and slightly less bound isovector 1 and isoscalar 1 h 
states. 

The pattern of relative binding energies between the channels may be understood from Fig. [5j The most deeply 
bound states occur in the isoscalar 1 h channel where the potential is repulsive near the origin but has a deep 
attractive well (due again to the isospin factor of 3) near lfm. The second most deeply bound states occur in the 
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r/fm r/fm 

FIG. 5: The potential, Eq. (|10[1 . plotted against r in the various 1" channels for h = 1 and A = lGeV. The dotted lines are 
the isovector potentials while the solid lines are the isoscalar potentials. The left panel shows the C = — potentials and the 
right panel gives the C = + potentials. 









Binding Energy / MeV 




State 


Isospin 


h = 0.8 


h = 0.9 


h = 1.0 


h = 1.1 


h = 1.2 


h = 1.3 


IS (0,1,2)— 





12 


20 


29 


60 


110 


160 


2S 




1.6 


3.6 


23 


39 


51 


65 


3S 






0.7 


6.7 


11 


15 


21 


IS (0,1,2)— 


1 


4.2 


8.8 


15 


21 


29 


38 


2S 








0.2 


0.7 


1.5 


2.8 


IS (0,1,2)"+ 





47 


67 


90 


120 


150 


180 


2S 




4.2 


8.1 


13 


19 


27 


35 


3S 




0.5 


1.6 


3.5 


6.1 


10 


14 


IS (0,1,2)"+ 


1 


0.1 


0.5 


1.6 


3.4 


5.9 


8.9 


2S 










0.1 


0.4 


0.9 



TABLE IV: Binding energies for D\D states with various isospins and charge conjugations; the binding energies are given for 
a few values of h in the range identified above and A = lGeV. 

isoscalar 1 channel where the • Tj term contributes a -3 factor and there is a deep attraction near the origin. We 
can see in Fig. [5] that the form-factor has reduced the magnitude of the first dip in the oscillating potential (around 
0.3 fm), making it smaller than the first bump (around 1.2 fm). This is why the isoscalar 1 ^ channel has deeper 
binding than the 1 channel. The isovector channels lose the isospin factor of 3, leading to significantly reduced 
binding in these channels. However, their relative binding is the same: the 1 potential retains the deeper attraction 
near 1 fm whereas the isovector 1 h channel loses the attraction around the origin due to form- factor effects. Hence 
the isovector 1 h channel is the least deeply bound when A = 1 GeV. Consequently the prediction of bound states 
in the isovector 1 h is the least robust. 

The situation is very different for the isovector 1 and isoscalar 1 h channels. In both of these channels the 
point particle potential is repulsive at the origin and they must bind in the first attractive well which is wlfm away 
from the origin. Therefore we expect these numerical results to be robust against finite size effects and a repulsive 
core. However, in the presence of intense regulation of the potential then, in this channel, the deep attraction being 
overwritten to strong repulsion with decreasing A, shown in Fig. [1] becomes a strong repulsion being overwritten 
to deep attraction. Therefore, we conclude that the existence of deep binding in these channels is a very robust 
result which should be insensitive to strong, short-range dynamics and totally independent of finite size effects of the 
potential, though both may contribute to deeper binding. 

The ranges of h and A which reproduce binding energies for the y(4260) and F(4360) (see Fig. @| are of particular 
interest. The case h=1.3 (Table ITVl) illustrates how it is possible to identify the IS and 2S 1"- respectively as Y~(4260) 
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and V(4360). In such a scenario it is possible that a third 1" state could occur around 4400 MeV. But of most 
interest is the prediction of a robust isoscalar exotic 1 h bound state in the vicinity of, or even below, the y(4260). 
If this exotic state is below the V(4260) then it may possibly be observed through V(4260) — > y(4200?) + 7. 

Table ITVl shows binding in both isovector 1~ channels. We therefore must reverse our previous concurrence @ with 
the conclusions of Ref. [ljjj: when subjected to a more complete analysis we find that a bound state may exist due to 
one pion exchange between D\D* in the isovector 1 _ channel. 

We find it interesting to note that the Z(4430) has a mass of 4433 ± 4 MeVfnj. Therefore if it were a D-JF 
molecule, it would have a binding energy of 4 ± 9 MeV. This binding energy is compatible with a charged partner 
of the 1 v isovector result for the range of h and A which reproduces the Y(4260) and V(4360). A more complete 
analysis than that provided here is necessary to make a definitive identification. However we find the possiblity that 
one pion exchange might provide a consistent description of the V(4260), V(4360), and the Z(4430) with physically 
reasonable parameters encouraging. 

In addition, we predict doubly charmed {D\D* as opposed to D\D*) isoscalar and isovector states degenerate with 
respectively the isoscalar and isovector D\D* states in C = — . We refer to Ref. j|| for a discussion of the signs 
involved. 



D. Bottom analogues 

In Table |V] we present the binding energies of some of the lowdying isoscalar B\B* states in relative S-wave 
with C = — , along with the analogous D\D states for comparison. Binding energies are given for a few values of 
h G [0.8,1.3] and A = 1 GcV. 







Binding Energy / MeV 




State 


h = 0.8 


ft = 0.9 


h = 1.0 


h = 1.1 


h = 1.2 


h = 1.3 


IS DjD* (0, 1,2) — 


12 


20 


29 


60 


110 


160 


2S 


1.6 


3.6 


23 


39 


51 


65 


3S 




0.7 


6.7 


11 


15 


21 


IS B{B* (0,1,2)— 


56 


93 


140 


190 


250 


320 


2S 


20 


29 


38 


49 


61 


75 


3S 


6.3 


9.8 


14 


19 


24 


30 



TABLE V: Binding energies for various isoscalar B\B and D\D states in L — with C = — ; the binding energies are given 
for a few values of h in the range identified above. 

The binding energies are generally deeper than in the charmed analogues. This is easily understood: the higher 
mass of the B mesons result in a lower kinetic energy. In general we predict analogous effects in the B analogs of the 
charmed system, subject to differences in the width which is experimentally undetermined for the B\. 

Similar effects may exist in the K system. However the phenomenology of the K\ is more complex and the heavy 
quark approximation is certainly inadequate. Together with the constraint on h implied by the width, this prevents 
us making quantitative conclusions, we only note the qualitative possibility that S-wave pion exchange may produce 
binding in the K system. 



V. DISCUSSION 

The results for binding energies, and even whether states bind at all, are sensitive to parameters, and also to more 
complicated (possibly more realistic) modelling of the strong interactions. 

We have focussed solely on the t-channel force from virtual pion exchange, specifically, the four fermion intermediate 
states in the Fock state. Therefore, we have taken only the real part of the potential and solved the Schrodinger 
equation and ignored the imaginary part arising from the exchange of a real, on-shell, pion. There are also s-channel 
forces arising from intermediate cc excited states. More immediately in our molecular approach there are intermediate 
states with a real pion, of form D*wD*. 

The ability of a virtual exchanged particle to be on-shell introduces an imaginary component to the matrix element 
and, hence, to the potential. The effect is to make the energy complex: the real part is taken as the binding energy 
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while the imaginary part is interpreted as the width of the state. That the on-shell intermediate state would manifest 
as a width seems natural as it represents a direct connection between the bound state and a possible decay channel. 

The picture is then that the D\ decays into a D*tt and the "would be" quasi- molecular bound state disintegrates, or 
even fails to form. Thus, we expect the on-shell pion contribution will endow any state produced by this mechanism 
with a width, or that it produces a non-resonant background which may obscure the signal. 

Within our approximations we find deeply bound meta-stable states. The D*irD* generates widths and background. 
Whether these states remain visible is then dependent upon the relative importance of neglected forces, such as mixing 
with cc or D*tvD*. In general it is difficult to calculate the impact of neglected effects, not least because strong 
interactions are complicated and we are approximating one particular force as dominant. 

If the y(4260) is an example of our states, then its visibility shows that Nature is kind, at least in the 1" channel. 
It has given a width of O(100)MeV and a visibility above background. It could be that this fortune is because a cc 
component drives the production, and the D\D* rearrangement then drives the ipirir signal. 

Thus the conclusion of this analysis is that while it is possible that a deep bound molecular spectroscopy with 
signals visible above a background can arise, it is not mandatory. However, as we have already noted, the appearance 
of V(4260) and Y(4360) are consistent with being the first two states observed in such a spectroscopy. The immediate 
test of this is to seek evidence for these states in DDtttttt. Unless there is some dynamical suppression, such channels 
must show strength if a D*D\ bound system is present. If this first test is passed, then a search for other transitions 
and evidence for analogous states in B* B\ would be warranted. In this latter case we note the apparent presence of 
anomalous state T(10.88)[l8| 

This and other phenomenological implications are the theme of the next section. 



VI. PHENOMENOLOGY 

We have studied the D*D\ molecules and found deeply bound states with I — 0, which are degenerate for the 
(0, 1, 2) channels. However, the number of potentially deeply bound states is very sensitive to parameters. Typically 
we anticipate the binding energies of the I — states to have the orders of magnitude as follows: IS O(10— 100)MeV; 
2S 0(1 - 10)MeV, with an exotic 1"+ between the IS and 2S levels. 

Further reasons to anticipate a rich spectroscopy are that this S-wave tt exchange also can occur for D — > D and 
the off-diagonal DD\ — > DqD* . The strengths for each of these in the heavy quark limit are identical. In practice 
there will be model dependent perturbations due to mass shifts and mixings; these are beyond the present paper and 
only merit study if the general features of our model show up in the data. 

In general: if S-wave pion exchange forms deeply bound charmed molecules comprised of D*D\\ D*Dq; DD±; DDq 
(and manifest charm analogs) , there will be a rich spectroscopy of states in the 3.9 — 4.5 GeV mass range. These 
can include states that are superficially charmonium, such as 7 = 0, h and 1 , as well as exotic J PC : and 
1 In addition there are also states with charmonium character but 1 = 1. States such as I = 0, h and 1~~ 
may contain cc in their Fock state and hence be produced at measurable rates; the other states have no such aid, 
but may be produced in radiative or strong transitions from higher lying molecules. Manifestly charm {D*D\, etc.) 
states are also expected and are degenerate with the C — — charmonium like states. The pattern and observability 
of these will depend on the detailed pattern of the spectroscopy. 

The states that are most amenable to experimental study are the 7 = 0, 1 . These occur in D±D*, and also can 
arise from the off-diagonal S-wave potential for DD\ — > DqD* . Hence there can be a rather rich spectroscopy in the 
1 sector. As there are apparently several states of varying statistical significance emerging in the data, we shall 
primarily focus on this channel here. 

The best established enigmatic structure in the 1 sector here is y(4260) which is seen in ipinr. Its typical hadronic 
width T(4260) ~ 90MeV implies that either ipirir is not the dominant decay channel or that 40 years of experience 
with the OZI rule and strong interactions is wrong. Given the nearness of the D(L = 0) + D(L = 1) thresholds which 
can be accessed in S-wave, rearrangement into ?/>7r7r at low momentum seems reasonable, and has been invoked as a 
qualitative explanation of these phenomena (19|. 

As D* and Dq — > Dtt, whereas D\ — > Dtttt, then if the dynamics are associated with the nearby DD\ and D*Dq 
thresholds, such as y(4260) being a DD\ molecule, or a hybrid cc that is dynamically attracted towards that 
threshold, then strength should be seen in the DDirir channels [19J, |20J. However, if the F(4260) is a D*D\ bound 
state, then the favored strong decay will be — > DD3tt in contrast to the aforementioned DD\ or D*Dq —> DD2tt. 

A preliminary report from Belle [24| sees no evidence for D*Dir in the F(4260) region. This disfavors DD\ and 
potentially also the D Dtttt channel. Thus by default, the possibility that the strength is driven by DDtttttt becomes 
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tantalizing. 

Thus an immediate consequence of this interpretation is that if the F(4260) is a D*D\ molecule, there must be 
significant coupling of the F(4260) — > DDtttttt that could exceed that to DDirir. More generally an unavoidable 
conclusion of this dynamics is that in the 1 sector the e + e~ — > DDtttttt channel has significant strength in 
the region of any D* D\ molecular states. Hence we urge measurement of the relative importance of the channels 
e + e~ — » DDtttttt and of e + e~ — > DDtttt (when, in the latter, D*D* has been removed). 

The depth of binding of the ground state with trial wave functions already suggested @ the tantalizing possibility 
that a radially excited state could also be bound. Numerical solutions of the Schrodinger equation confirmed that 
this is likely to be the case in the range of models discussed here. 

The excitation energy for radial excitation of a compact QCD cc state is 0(5OO)MeV; it takes less energy O(100)MeV 
to excite the extended molecular system which has no linearly rising potential. The spatial extent of the molecular 
2S system is significantly greater than that of cc hadrons. The rearrangement of constituents leading to final states 
of the form if; + light mesons then rather naturally suggests that the lower (radial) states convert to ipTnr ( %P'tttt) 
respectively. In this context it is intriguin g th at there are states observed with energies and final states that appear to 
be consistent with this: F(4260) — > ipTnr [2lj j and the possible higher state Y(4360) — > ^'tttt [H, [23| are respectively 
170MeV and 70 MeV below the D*(2010)Z>T(2420) combined masses of 4430MeV. Here again, for a D*D~i molecular 
state, we would expect significant coupling to DDtttttt. 

If these states were to be established as members of molecular systems, one could tune the model accordingly. 
Further, this could be an interesting signal for a D*D\ quasi- molecular spectroscopy with transitions among states 
that could be revealed in, for example, e + e~ — > tpjjTnr. Indeed, if we identify 1S(4260) and 2S(4360), then we expect 
the exotic 1 h to occur in the vicinity of the F(4260). Given that lattice QCD finds activity for a hybrid cc* signal in 
this channel in this region, one should now actively search for evidence. A clear signature is that the 1 ^ hybrid will 
couple to DDtttt in either the D*Dq or DD\ combinations; looking for the presence of strength in e + e~ — !• ^DDtttt 
which does not include D*D* should thus be a primary endeavor. The absence of such a channel could have far 
reaching implications for theory. 

While our discussion has centered on charmonium, the remarks hold more generally. Since the attraction of the 
potential depends only on the quantum numbers of the light qq, it follows immediately that the flavor of the heavy 
quarks is irrelevant, at least qualitatively. Hence we expect similar effects to occur in the bb and ss sectors. It has 
been noted that T(10.86)GeV appears to have an anomalous affinity for T7T7T 25j. This state is ~ 130 MeV below 
B*B\ threshold. In the cftmr channel there is an enhancement at 2175MeV[26j]. This is approximately 125MeV below 
the K* Ki(lA00) threshold. This is consistent with the K*K\ spectroscopy; however, as commented earlier, analysis 
here is less reliable, as the heavy quark approximation fails, and the phenomenology of the JsTi(1270; 1400) pair is 
more complicated [271 l28j . 

The primary test for this picture is that if the states in the 4 to 4.5 GeV region are deeply bound D*D\ spectroscopy, 
then their decays in charm pairs must show strength in the DDtttttt channels. The energy dependence of this channel 
and that of DDtttt (with no D*D*) can reveal the mixings between D*D\ and DD\/ D* Dq molecular systems. The 
presence of exotic 1 '"is also expected. 



VII. CONCLUSIONS 



In general we find that deeply bound molecules in the D\D* system should occur as a result of tt exchange in S-wave, 
leading to a potentially rich spectroscopy. Whether such states are narrow enough to show up above background is a 
question that experiment may resolve. We note however that the emerging data on the 1 states known as F(4260) 
and F(4360) are consistent with being examples of these molecular states. The immediate test is to verify if the 
prominent channels with manifest charm in this mass region are DDtttttt. If this is confirmed, then more detailed 
studies will be merited, in particular searches for an exotic 1 h in the vicinity of 4.2GeV. This state could be produced 
via V(4260) — > (1 h ) + 7, and/or be revealed in 1~ + — > ipj. 

Table HVl with h=1.3 shows a possible spectroscopy consistent with the Y(4260) and Y(4360) as the IS and 2S 1~~ 
states. In this case, the exotic states expected are 1=0 0~~ also at 4260 and 4360 (in both charmoniumdike and 
manifestly charm channels); the isoscalar 1 h at 4250 and 4395; and also 1=1 "charmonium" states, including 1" 
at 4390. 

As long as one picks and chooses which datum one will fit, it is possible to fit it in a molecular model. A reason is 
that binding energies are very sensitive to parameters that are not well determined elsewhere. Thus a model designed 
to fit a single state has limited appeal. The more relevant test is whether a group of states share a common heritage, 
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and their production or decay properties reveal the underlying molecular structure. In the particular case here, one 
can fit the masses and decay widths in tetraquark, hybrid and molecular models. As such the existence of these states 
does not discriminate among them. 

However, the pattern of J PC and the decay channels differ. Thus the sharpest tests of their dynamical structure 
appears to be in the decay branching ratios. Hence, for example, the K(4260) as a cscs tetraquark would couple to 
D S D S ; a hybrid or molecule associated with DD\ threshold would be expected to appear in DD\ — > DDirir; molecules 
associated with the D*D\ threshold by contrast would have significant strength in the DDiririr channels. Thus the 
decay branching ratios of states seem likely to be sharper indicators of their dynamical nature than simply their 
masses. 

If our hypothesis is correct, we expect significant strength in the e + e~ — > DDiririr channels in the 4 to 5 GeV 
region. Such evidence may already exist among the data sets for e + e _ annihilation involving ISR at BaBar and Belle. 
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